ASYMPTOTICS OF CHARACTERISTIC POLYNOMIALS 
OF WIGNER MATRICES AT THE EDGE OF THE SPECTRUM 



H. ROSTERS 

Abstract. We investigate the asymptotic behaviour of the second-order corre- 
lation function of the characteristic polynomial of a Hermitian Wigner matrix at 
the edge of the spectrum. We show that the suitably rescaled second-order corre- 
lation function is asymptotically given by the Airy kernel, thereby generalizing 
the well-known result for the Gaussian Unitary Ensemble (GUE). Moreover, 
we obtain similar results for real-symmetric Wigner matrices. 



1. Introduction and Statement of the Main Results 
Let Q be a fixed probability distribution on the real line such that 

JxQ(dx)=0, a := J x 2 Q(dx) = 1/2, b := J x 4 Q(dx) < oo , (1.1) 

and for any N = 1,2,3,..., let := (Xij )ij=i,...,jv denote the associated 
Hermitian Wigner matrix of size N. This means that 



'X^+iXP? fori<j, 
V2X£ e for i = j, 

X]\ e -iX]f for i>j, 



where {X^ c | i < j} U {X}f \i < j} is a collection of i.i.d. real random variables 
with distribution Q. The second-order correlation function of the characteristic 
polynomial of the random matrix Xn is defined by 

f N (n,v):=E(D N (n)D N (v)) (p,u€R), 

where Dat(A) := det(Xj^ — A). We are interested in the asymptotic behaviour of 
/aK/UVj^aO as ./V — > oo, for certain sequences (/uv), (vn) which will be specified 
below. Furthermore, the correlation coefficient of the characteristic polynomial of 
the random matrix Xn is defined by 

. E((gg00 - EggQQ) {Dn{u)-ED n (u))) 
^E(D N (p) - ED N (fi)Y ^E(D N (u) - ED N (v)Y 
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In the special case where Q is the Gaussian distribution with mean and 
variance ^, the distribution of the random matrix Xjy is the so-called Gaussian 
Unitary Ensemble (GUE) (see e.g. Forrester [Fb] or Mehta |Me| but note that 
we work with a different variance). In this case, it is well-known that 

f N (H, u) = V^N\ e (M 2 +- 2 )/4 K N+ i(fi, v) , (1.2) 

where 

ti v^(fc-l)! 

and the p^ are the monic orthogonal polynomials with respect to the weight 
function e~ x I 2 (see e.g. Chapter 4.1 in Forrester [E>]). Thus, up to scaling, 
the pk coincide with the Hermite polynomials (as defined in Szego [Sz]), and it is 
possible to derive the asymptotics of the second-order correlation function /jv from 
the well-known asymptotics of the Hermite polynomials (see e.g. Theorem 8.22.9 
in Szego [Hi]). More precisely, one obtains the following (well-known) results 
(see also Chapter 4.2 in Forrester [Fbj): For £ G (-2, +2) and any /j,, v e R, 

Jim c' N f N (v^V£ + fi/VN (O, VNZ + v/VNq(0) = S(/x, u) , (1.3) 

iv— >oo \ / 

where c' N := (V^NIN 1 / 2 e (0 eMlm 2 +U^+^/ q(0))~\ e (0 ■= 
and 

sin-7r(u — u) , 
S(H,U):= ^ / ■ (1.4) 

For £ = +2 and any ^i/el, 

lim 4 /tv f 2\ZjV + /i/iV 1/6 , 2ViV + i^/iV 1 / 6 ) = A(p, u) , (1.5) 

where := N\ N 1 ^ exp(2iV + (jj, + ^iV 1 / 3 ))- 1 , 

= Ai M Ai M -Ai M Ai Mj 
/x — V 

and Ai denotes the Airy function (see e.g. Abramowitz and Stegun |AS| ). 
By symmetry, a similar result holds for £ = — 2. The functions in (jl.4p and (|1.6p 
are also called the sine kernel and the Airy kernel, respectively. Furthermore, 
it is well-known that the eigenvalues of a random matrix Xn from the GUE are 
distributed roughly over the interval [— 2y/~N, +2y/~N]. That is why the results 
(II. 3p and (|1.5p are also said to refer to the bulk and the edge of the spectrum, 
respectively. 

Recently, Gotze and Kosters |GK| have shown that the result (|1.3p for the bulk 
is (almost) "universal" in the sense that it holds not only for the GUE, but also 
(with minor modifications) for more general Hermitian Wigner matrices as intro- 
duced at the beginning of this section. More precisely, under the assumption (jl.ip . 
we have 

lim 4 f N (ViV£ + li/VNg®, y/N£ + »/VNq(0) = exp(6-f ) S(/i, v) (1.7) 

for all £ £ (—2, +2) and all fi, v G R, where Cjy, g(£) and S(/i,z/) are the same as 
in (TOD. 
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Therefore, it seems natural to ask whether the result (|1.5p for the edge can also 
be generalized to more general Hermitian Wigner matrices. The main purpose of 
this paper is to answer this question in the affirmative. More precisely, our first 
result is as follows: 

Theorem 1.1. Under M.l}) . we have 

lim c" N f N (2VN + fi/N^^^VN + =exp(6-|)A(//,i/) (1.8) 

N— >co V / 

for all v G R, where c" N and A(/u, v) are the same as in i f 1.5\) . 
Corollary 1.2. Under we have 

lim a N foy/N + /m/N 1 /®, 2^N + u/N 1 /®) = A^ j == (Lg) 



v / A(^,/i) a/A(^ i') 



/or a// /u, z/ G 



Moreover, it turns out that similar results hold for real-symmetric Wigner matrices. 
Let Q be a fixed probability distribution on the real line such that 

/,%*,) = <>, S :=/^) = l, i^/^Q^Xoo, (1.10) 

and for any N = 1,2,3,..., let Xn '■= (-Xy)ij=i,...,jv denote the associated 
real-symmetric Wigner matrix of size N. This means that 

'Xff for i<j, 

V2X% C for i = j, 
y X^ for i>j, 

where {Xf- e \ i < j} is a collection of i.i.d. real random variables with distribution 

Q. Then, similarly as above, the second-order correlation function and the cor- 
relation coefficient of the characteristic polynomial of the random matrix Xn are 
defined by 

f N (fi,u) :=E(D N (p)D N (y)) (/i, v eR) 

and 

~ , , E((D N (ri-ED N (ri) (D N (v)-V,D N (v))) 

a N {ji,v):= = (/i,i/GM), 

y/E(D N {fi) - E-Djv(aO) \]HD n {u) - ED N (v)) 

respectively, where now Dn(A) := det(X/v — A). 

Following the approach by Gotze and Kosters [GKj . Kosters [Ko] recently 
showed that under the assumption (jl.lOp . we have 

lim d' N f N ( v^e + ///VJVeCO. ^ + v/"/NQ(0) = exp(%^) T(fi, u) (1.11) 



for all £ G (—2, +2) and all G R, where £>(£) is the same as in (jl.3p . 
dJv := (v^FiV!iV 3/2 «(0 3 exp^iV^ + KM + ^CMO)) -1 , and 

2sin7r(^- v) 2costt(^ - t/) 

T(A*,v := 7 7 v» — • (1-12) 
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Thus, we also have universality (in the same sense as above) in the bulk of real- 
symmetric Wigner matrices. In the special case where Q is the Gaussian distri- 
bution with mean and variance 1, the distribution of the random matrix Xn is 
the so-called Gaussian Orthogonal Ensemble (GOE) (see e.g. Forrester |Fo] or 
Mehta [Me] but, again, note that we work with a different variance). In this case, 
([blip had been obtained previously by Brezin and Hikami [BH2]. 

Our second result shows that the result (jl.lip admits an analogue for the edge 
of the spectrum, too: 

Theorem 1.3. Under il.lO\) , we have 

lim d% J N (2VN + /i/iV 1/6 , 2VN + u/N 1 ^) = exp(%^) M(p, v) (1.13) 

for all {j,,veR, where d" N := N\ N 1 / 2 exp(2iV + (/x + z^JV 1 / 3 ))- 1 , and 

(/z + I /)Ai( M )Ai(i/)-2Ai'(/z)Ai>) 2Ai( / u)Ai>)-2AiV)AuV) 
U{fi, v) := -2 V 



(1.14) 



Corollary 1.4. Under il.l(J\) . we have 

lim a N (2VN + v/N^^VN + v/N 1 ^) = B( ^ } - (1.15) 

for all fj,, v E M. 

In the special case of the GOE, 13j) can already be found in Brezin and Hikami 
[BH3] . 

It seems interesting to note that the functions § and T arising for the bulk of 
the spectrum are related by the identity 

n^y) = (-^—.(i-.-i:))H^y) 



x — y \dy dx J 

and that the functions A and B arising for the edge of the spectrum are related by 
the analogous identity 

see also Brezin and Hikami [BH3] . (To check the latter identity, use the fact that 
the Airy function Ai(z) satisfies the differential equation Ai"(z) = zAi(z).) 

Also, observe that in all the cases previously mentioned, the precise choice of 
the underlying distribution Q or Q enters into the asymptotic behaviour of the 
second-order correlation function of the characteristic polynomial only as a multi- 
plicative factor depending on the fourth cumulant. Thus, the results are essentially 
"universal" . 

Let us mention some related results from the literature. It is well-known that at 
the edge of the spectrum of Wigner matrices, we have universality also for the cor- 
relation function of the eigenvalues themselves (see Soshnikov [So]). In contrast 
to that, for the bulk of the spectrum of Wigner matrices, only partial results are 
available in this direction (see JOHANSSON [Jo] ) . Furthermore, in the special cases 
of the GUE and the GOE, several authors have investigated the averages of more 
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general products and even ratios of characteristic polynomials (see e.g. Brezin 
and Hikami [BHTl lBH2| IBH3] . Fyodorov and Strahov [FS], Baik, Deift 
and Strahov [BDS], Strahov and Fyodorov [SF], Akemann and Fyodorov 
|AF| . Vanlessen [Va], Borodin and Strahov |BS] ). Even more, at least in 
the Hermitian setting, some of these results have been shown to be "universal" 
in that they continue to hold (with some modifications) for the class of unitary- 
invariant ensembles. For Wigner matrices, however, less seems to be known in 
this respect. 

Acknowledgement: I would like to thank Alexander Soshnikov for the suggestion 
to study the problem at the edge of the spectrum. 



2. Outline of the Proofs 

To prove Theorems 1 1.1 1 and 1 1.31 we will start from the fact that for fixed 
the exponential generating functions of the sequences fpf(pb,v) and /jv(/^) v) are 
given explicitly by 

^ X N exp(^- T ^-i(/i 2 + ^ 2 )- T ^ + (6-f)x 2 ) 

t/ N ^ V) M = (1-S)3/2.(1 + *)!/* (N < ^ 

(see Lemma 2.3 in Gotze and Kosters |GK| ) and 

~ X N exp(^. T ^-i(^ + ^). T ^ + (^)x 2 ) 

tJ NM m = (i-*)s/ 2 .(i + *)V 2 ^ < V 

(see Lemma 2.3 in Kosters |Koj ). respectively. This fact opens up the possibility 
to study the asymptotic behaviour of the second-order correlation functions by 
evaluating appropriate contour integrals of their exponential generating functions. 
In fact, this strategy was already used by Gotze and Kosters |GK] and Kosters 
[K6] to obtain the above-mentioned results for the bulk of the spectrum. Here we 
carry out a similar analysis for the edge of the spectrum. 

Since it does not require any additional efforts, it seems convenient to evaluate 
the values 



(2.1) 



Nl ' 2ni 7 7 (1 - z )«+(i/2) . (i + z )i/2 z n+i 

for arbitrary a > and 6* 6 K, where 7 denotes a contour around the origin. By 
the foregoing, the case a = 1 corresponds to the Hermitian case and the case a = 2 
corresponds to the real-symmetric case. 

We remark in passing that for a general parameter a > 0, it is not hard to see 
that the exponential generating function under consideration can be interpreted 
as that of the second-order correlation function of the characteristic polynomial of 
a random matrix from (a rescaled version of) the tridiagonal beta ensemble intro- 
duced by Dumitriu and Edelman [DE] with a = 2/f3. For this interpretation, 
one should set b* := 0. 
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For a > and fj,, v G R, put 
(Q) _ J_ f+°° exp - |fr + - |Qx - u) 2 /(l-iu)) 

^' ^ • 47r 3/2 y_ oo ( X _ iu )a+(l/2) U ■ 

(2.2) 

In the next section, we will prove the following results: 
Proposition 2.1. For any a > 0, b* G R and (j^eB, u/e /iai>e 

lim (V^iV! AT^- 1 )/ 6 exp(2iV + (// + i/)^ 1 / 3 ) 

• /^ ) (2v / iV + /i^ 1/6 ,2v / iV + ^" 1/6 ) =exp(6*)/(°)(/i,i/). 
Proposition 2.2. For any a G N and /i^Gl, we /iaue 

^^) = (^(^-|;)) w K)ak.: 

where (-)^ denotes the a-fold application of the given differential operator. 
(For fj, = v, consider the continuous extension of the right-hand side.) 

It is straightforward to check that 



and 



^(«.») - (!" i)) <2> (AiWAifo)) -B(,,„). 

(For the second identity, use the fact that the Airy function Ai(z) satisfies the 
differential equation Ai (z) = zAi(z).) Thus, in view of the preceding comments, 
it is immediate that Theorems 11.11 and 11.31 which correspond to the special cases 
a = 1 and a = 2, follow from Propositions I2.l l and 12.21 
To prove Corollaries 11.21 and 11.41 observe that 

, v f N (ji,v)-®D N (ji)ED N (v) 

yJfNfan) - $Dn{h)) yJfNiy,v) - (^D N (u)Y 

in the Hermitian case and 

f N (n,v)-ED N (ji) ED N (u) 



y/ JnM - {ED N (v)) 2 y/ >jv (!/,!/) - {ED N (u)) 2 

in the real-symmetric case. Moreover, it is not difficult to see that 

ED N (X) = ED N (X) = g N (X) := (-1) N 2~ N / 2 H N (X/V2) (A G R) , 

where Hn(x) is the A^th Hermite polynomial (see e.g. Section 5.5 in Szego [Sz]). 
Thus, in both cases, the expectation of the characteristic polynomial is given 
by the same function <7jv(A). Therefore, to deduce Corollaries 11.21 and 11.41 from 
Theorems 11.11 and 11.31 respectively, it will be sufficient to show that gN{ft) 9n{v) 
is asymptotically negligible in comparison to /at(/U, v) and /j\r(/x, v). 
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Slightly more generally, we will consider the case of an arbitrary parameter a > 
and investigate the asymptotic behaviour of 

°N (M. V) ■= I r— ( 2 - 3 ) 

yfu i^v) -9n{v) 2 yfu iy, v)- gN (v) 2 

for any fj,, v E R. In the next section, we will show that Proposition 12.11 entails 
the following result: 

Proposition 2.3. For any a > 0, b* G R and fj,,u el sitc/i i/iai /( Q )(/i,/i) > 0, 
/( a ) (i/, z/) > 0, we /iawe 

Jim a^\2VN + ^N-y\2^ + vN-y & ) - l(a) ^ v ) 



Furthermore, we will prove the following: 

Proposition 2.4. For any a G N, we have I^ a \x,x) > for all i£R. 

In view of the preceding comments, it is obvious that Corollaries 11.21 and 11.41 
which correspond to the special cases a = 1 and a = 2, follow from Propositions 
EH and El 

We remark in passing that for a general parameter a > 0, crjy (/^j z/ ) can be inter- 
preted as the correlation coefficient of the characteristic polynomial of a random 
matrix from (the rescaled version of) the tridiagonal beta ensemble (with a = 2/(3), 
since in this setting, the average of the characteristic polynomial is also given by 
the Hermite polynomial (see Theorem 4.1 in Dumitriu and Edelman [DE] ) . 



3. The Proofs 
Proof of Proposition \2.1\ Fix fj,, v £ R. We have to evaluate 

A a Vjv,i/jv) j_ r ex P (jwn • ijg - j04 + v%) • izg + 

AH ~ 2vri ./ 7 (1 - z) Q +( 1 /2) . (i + z )i/2 ' 

(3.1) 

where ■= 2N 1/2 + fiN~ 1/G , v N := 2N 1/2 + vN~ 1/6 , and 7 denotes a contour 
around the origin (which will be chosen further below). 

Setting £jv := (fJ>N + ^at)/2 and t/at := (//jv — ^iv)/2, (13. 1|) may be written as 



1 /■ ex P (~2?JV • T+i ~Vn • — + & 



dz 



exp (le 2 , + ^) ■ — ^ {l _ z)aH m . (i + z y/2 iivn ( 3 - 2 ) 

with 

? N = AN + 2{ii + v)N 1 '* + \( l i + v) 2 N- 1 ^ and j& = - ufN' 1 ^ . 



In particular, the leading exponential factor in (|3.2|) is asymptotically the same 
as that in Proposition 12.11 
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Thus, to complete the proof of Proposition 12.11 it suffices to show that 



N -{2a-l)/& r exp v ^ 

N^oo (27r) 3 / 2 i 7 7 



If 2 . ]-£ _ 2 . JL -L. frV 2 

SJV 1 + 2 ^TV 1-2 + z 



(l_ z )a+(l/2) .(l + Z )l/2 



exp(&*) 

4^ 3 /2 



+oo 



exp (/too (1 - 

(1 - iu)«+(V2) 



where 



hoo(z) := -kz 3 -\(n + v)z-\(n- vfjz . 



du , (3.3) 



(3.4) 



Similarly as in the proof of the main theorem in [GK] . the basic idea is that 
the main contribution to the integral in (13. 3p comes from a small neighborhood of 
the point z = 1. Let 

1 



1-z 



h N {z) : =-ffi.l— ^-rf N .-^- + b*z 2 -{a + \)\og{l-z)-\\og{l + z)-N\ogz 

(where log denotes the principal branch of the logarithm) and 

7iv(t) := (1 - iV~ 1/3 ) exp (it) , -ir < t < +vr . 
Then the left-hand side in (13.31) can be rewritten as 



N -(2*-l)/6 [ exp( y -±e N -TTz- r l 2 N-Th + b * z 
(2vr) 3 / 2 i 



7jv 



(1 - Z )«+(l/2) • (1 + z)V2 



dz 



( 27r) 3/2 z - / «p(M*))- 



jy-(2a-l)/6 f+Tr 



7JV 



(2tt) 3 / 2 



exp(/iAr(7Ar(t))) (it. 



Put/jv,i(a) := (-a^" 1 / 3 ;+aiV- 1 / 3 )and/7V,2(a) := (-vr, +vr)\(-aiV- 1 / 3 ; +aiV- 1 / 3 ), 
where a > 0. We shall prove the following: 

Claim 1: For any fixed a > 0, 

jy-(2a-l)/6 



lirn 

N^oo (2vr) 3 / 2 



, / + ^ m. ex P 6 ) / exp (/too (1 - iu)) 
exp(h N 7Ar (t ))) c/t = , / — . , a+(1/2) du . 

I N1 (a) 47T 3 / 2 J_ a (l-m)«+lVZJ 



Claim 2: For any 5 > 0, there exists a constant do > such that for all a > ao, 

Ar -(2a-l)/6 



for all iV G N sufficiently large. 



exp(/ijv(7iv(t))) dt 



< 5 



Before turning to the proofs, let us show that Claims 1 and 2 imply (|3.3 
Observe that 

30 exp (/too (1 - i«)) 



(1 - iu)«+( 1 /2) 



du < oo . 
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Thus, for a > sufficiently large, we have not only the conclusion of Claim 2, but 
also the inequality 



exp(&* 



47r 3/2 



exp(/t oc (l - iu)) 

-a,+a) (l-i«)«+(V2) 



< 5. 



Hence, in combination with Claim 1, it follows that 

N -(2a-l)/6 f +* exp(6*) 

exp(/i 7V ( 77V (t))) dt 



< 



(2tt) 3 / 2 

Ar -(2a-l)/6 



4^3/2 



+oo 



exp(/ioo(l - iu)) 
(1 - m)«+(V2) 



(2vr)3/2 



%,2(l) 



exp(/iAr(7Ar(t))) 



+ 



+ 



N -(2a-l)/6 



(2tt) 3 / 2 
exp(6* 



exp(6*) /" +a exp(/ioo(l - iu)) 



exp(/i JV (7j V (t))) dt . 3/2 / . , .„ 

I N1 (a) ^ ' J -a (1 - lu) a+ WV 



du 



exp(/ioo(l - zit)) 



du 



47T3/2 7 RX( _ 0)+0) (l_i«)«+(V2) 

for all N £ N sufficiently large. Since 5 > is arbitrary, this proves (|3.3 
Proof of Claim 1: First of all, substituting t = uN^ 1 / 3 , we have 



N 



-(2a-l)/6 r+aN-V 3 



jy-(2a+l)/6 r+a 



exv(h N (j N (t))) dt - 
oJV-i/3 (27r) 3 / 2 



exp 



(/i JV ( 7i v(«^- 1/3 )))d«. 



(2vr) 3 / 2 

We will determine the asymptotics of /iat(7at(u./V _1 / 3 )) as iV — > oo, for u G [—a, +a]. 
The O-bounds occurring in the sequel hold uniformly in this region. To begin with, 
note that 

z := 7 Jv(niV- 1 / 3 ) = (1 - N ^yuN-^ 

= (1 - iV~ 1/3 )(l + iuN~ 1/3 - \u 2 N- 2/z - \ivpN~ 1 + 0{N~ 4/3 )) 

= !_(!_ m)JV-V3 _ (iu + l M 2 )Ar -2/3 + ( l u 2 _ l iu 3 )JV -l + 0(^-4/3) 

Therefore, for iV sufficiently large, we have the following approximations: 



2 ^ T+z 

§&(j(l " *) + 1(1 " + |(1 " zf + 0((1 - *) 
-1 (4iV + 2(/i + z^)7V 1/3 + 0(A^~ 1/3 )) 

• (1(1 - iu)N-^ 3 + iAf- 2 / 3 + (1 + \iu + |n 2 - ^to 3 )^" 1 + 0(N^ 3 ) 
-(1 - m)iV 2 / 3 - \N 1 ' 3 + (i(l - iuf - !(// + - in) - l) + ©(A^ 1 / 3 ) 



1 - z 



-Vn 



1 - z 



(^-^- 1/3 )-(T3i^vr( 



1 + 0(iV 



-l/3> 



-i( M -^) 2 /(l-^).(l + 0(iV-i/3)) 
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b*z 2 = b* + 0{N~ 1 / 3 ) 
-(a + \) log(l -z) = +\{a + \) log N-{a + \) log(l - in) + 0(iV- 1/3 ) 
-±log(l + z) = -ilog2 + 0(iV- 1 / 3 ) 

-iVlogz = -N (iuN- 1/z + log(l - iV~ 1/3 )) 

= -iv (mj\rV3 - iv- 1 / 3 - ±Ar 2 / 3 - £ jv 1 + o(iv- 4 / 3 )) 

= (1 - m)iV 2 / 3 + iiV 1 / 3 + i + 0(iV- 1/3 ) 

Putting these approximations together, the terms of highest order cancel out, and 
we end up with 

IiN^NiuN- 1 / 3 )) = i (q + \) log TV + (^(1 - iu) 3 - + u){l - iu) 

- W - - *«) - (<* + ±) log(l - i«) + 6* - i log 2) + ©(iV- 1 / 3 ) . 

Hence, by an application of the dominated convergence theorem, it follows that 

jy-(2a+l)/6 r+a 



lim 



/ a e X p(/ lJV ( 7J v(niV- 1 / 3 ))) du 

J — a 



iV^oo (2tt) 3 /2 

exp(6*) /- +a ex P - + iu ) ~ W ~ ^) 2 /(i - »«)) 



-a 



4vr 3 /2 y_ a (l-m) Q +(V2) 

and Claim 1 is proved. 

Proof of Claim 2: By symmetry, it suffices to consider the interval (aiV -1 / 3 , ir). 
Write the integral as 



f'TV 

JaN- 1 / 



(2vr)3/2 J aN - 1/3 (l _ 7JV ( t ))«+(i/2) . (1 + 7JV (t))i/2 7JV (f)iV • 

Since £/v and i]n are real and | 7 jv(t)| < 1, this is clearly bounded by 
AT-(2a-D/6 r exp . Re (±zrf) _ ^ . Re + ^ 



(2vr)3/2 7 aAr _ 1/3 |l-7Jv(t)| a+ ( 1 /2).|l + 7Ar (t)|V2 |7iv(t)r' 

In the following, let K, Ki,K 2 > denote constants which depend only on a, b*,fj,,v 
(which are regarded as fixed) but which may change from occurrence to occurrence. 
Then we have 

S%>4JV- KN 1/3 

and 

( l-lN(t) \ _ 2N-^-N-^ 
6 Vl + 77v(i)/ (2 + 2cost)(l-iV- 1 /3) + jV-2/3 
(as follows from a straightforward calculation) and therefore 



1*2 
2 C 



1 - 7iv(i) \ 4iV 2 / 3 - 2N 1 / 3 - K 



1 + 7J\K*)/ ~ (2 + 2cost)(l- AT-V3) + jV-2/3 
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Since 



| 7J v(t)r = exp(iVlog(l-iV- 1 / 3 )) 

> exp(N(-N' 1/3 - ijV-2/3 _ KN- 1 )) 
= exp(-iV 2 / 3 - \N 1 / 3 _ K) , 



it follows that 

exp ^ 2 ^v Ke ^i+7jvW 



|7(*)I 



N 



( AN 2 ' 3 - 2JVV3 _ K x Ar2/3 lArl/3 ' 

^ " (2 + 2co S t)(l-iV-V3 ) + j v- 2 /3 + + ^ 1/3 + * 



< exp 



(2-2 cos i)A^ 2 / 3 - (1 - cos tjiV 1 / 3 - IT* % 
(2 + 2 cos i)(l - iV~V3) + AT-2/3 



for some constant -ftT* > 0, say. Now let e > denote a constant such that e 2 t 2 < 
1 - cost < \t 2 for all i G [0,vr]. Then, for a > e~ l V2K*, N 1 / 3 > maxje" 2 , avr" 1 } 
and t 6 [aiV _1//3 , tt], we have 

(2 - 2cost)iV 2/3 - (1 - cosi)iV 1/3 - K* 

> ^2^2/3 + { l £ 2 t 2 N 2/3 _ 1^1/3) + (1^2^2/3 _ ^ > ^2/3^2 

and therefore 



exp ^ 2 ?AT ±te ^l+7ivW 



l7(*)r 

< ^p f-7 w iy2/W / „ ^ < exp ( -N 2 / 3 e 2 t 2 /A) . 

~ \ (2 + 2cosi)(l-iV- 1 /3) + jv- 2 / 3 y / ~ V / 

Thus, since |1 ±7jv(i)| > iV -1 / 3 , the integral under consideration is bounded by 

O-M/e f exp(-AT 2 / 3 g 2 t 2 /4) 

AiV-V3 |l- 7 iv(t)| a+1 / 2 -|l+7^(i)| 1/2 ' 
<^AT-(2-i)/6 r eM-N^m 

- y aJV -va iv-(2 Q+ i)/6 

= KN 1/3 exp(-N 2/3 e 2 t 2 /A) dt 

roo 

<K exp(-e V/4) du . 

j a 

Obviously, this upper bound can be made arbitrarily small by picking a and N 
large enough. This proves Claim 2. 

The proof of Proposition 12.11 is complete now. □ 
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The basic ingredient for the proof of Proposition 12 , 21 will be the following integral 
representation for the product of two Airy functions: 

Lemma 3.1. For any x,y £ C ; 

w U . n 1 f exp (±z 3 -i( x + y)z- \{x - y) 2 / z) 

Ai(x) Ai(y) = ^ dz , (3.5) 

where C denotes some (unbounded) contour from ooe -7 ™/ 3 to ooe +7 ™/ 3 . 

In the special case where y = ±x, this result can already be found in Reid |Re| . 
For the convenience of the reader, we give a detailed proof of Lemma 13.11 



Proof. We start from the following well-known integral representation of the Airy 
function: 



Ai(x) = / expfiz 3 — xz) dz . 

2m J c 



A standard application of Cauchy's theorem shows that the contour C can be 
deformed into the contour 1 1— > 1 + it, t £ R. Thus, we obtain 

Ai(x) = — / exp(iz 3 - xz) dz . (3.6) 
2m Ji-ioo 

Observe that the resulting integral exists in the Lebesgue sense, since we have 

|exp(i(l + itf - x(l + it))\ < exp(± - t 2 + |x|(l + \t\)) 

for any t 6 E. 

It follows from (13.61) that 



Ai(x)Ai(y) = J J exp(|(l + itf - x(l + it)) exp(|(l + iuf - y(l + iu)) du dt . 

Substituting (t,u) = (^(v + w), \{v — w)) and doing a small calculation, we find 
that 

Ai(x) Ai(y) = -^JJ exp(±(2+ivf-l(2+iv)w 2 -±{x+y)(2+iv)-±(x-y)iw) dw dv . 
Using the well-known relation 

/+oo rz 
exp(-aw 2 - bw) dw = ^L e b 2 /^ 
-oo V a 

(where Re (a) > 0) , it follows that 

1 f exp(^(2 + ivf - i(x + y){2 + iv) - \{x - y) 2 /{2 + iv)) 



or 



AV w V , 1 [ 2 + io ° exp(^z 3 - |(a + y)z - \(x - y) 2 / z) 
Ai(x) Ai(y) = — ^ / dv . 



^' 2 i y 2 -ioo z^ 2 

By another application of Cauchy's theorem, the contour v i— ► 2 + iv, v £ R, may 
be deformed back into the contour £. □ 
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Proof of Proposition \2.S[ Replacing the contour £ in Lemma 13.11 by the contour 
t i— ► 1 + it and substituting t = —u, we have 

1 f + °° exp - \(x + y)(l-iu) - \{x - y) 2 /{l-iu)) 
Al (x)A l(y) = i -^y_ oo ^— du. 

By means of abbreviation, write E(x, y, u) for the numerator inside the integral. 
Then, for any a > 0, we have 

d ( f E(x,y,u) 



du 



dy \J (1 - iu) a J J ( 
Off E(x,y,u) \ f E(x,y,u 



a-tt (-K 1 -™) + sO c -y)/( 1 -*0) du - 



dx \J (1 — iu) c 
and therefore 



= / H(l-<«)-i(x-y)/(l-it*)) du, 



1 ^ <9 d\\ f f E(x,y,u) d \ _ f E(x,y,u) 



x — y \dy dx J J \J (1 — ui) Q / J (1 — iu) a+1 
The assertion of Proposition 12.21 now follows by induction. □ 

Proof of Proposition^^ Fix /j,,u £ R, and put /j^ := 2iV 1 / 2 + fiN' 1 / 6 , 
:= 2N 1 / 2 + vN^ 1 ^ . Using well-known results about the asymptotic properties 
of the Hermite polynomials (see e.g. Theorem 8.22.9 (c) in Szego |Szj). we find that 
the function g N {\) given by g N {\) := (-1) N 2~ N / 2 H N (X/V2) (A £ R) satisfies 

= 2~ N H N (V2N + fiN- 1 / 6 /^) H N (V2N + vN' 1 ^ j y/2) 

= O (2- N exp(AT + M ArV3) 2 w/2 ^,1/2 ^-1/12 exp(Ar + ^1/3) 2 N/2 N] i/2 N -l/l2^ 
= O (n\ N- 1 / 6 exp(2N + + zv)N 1/3 )' 



= o (N\ Af^- 1 )/ 6 exp(2AT + (fi+ ^N 1 / 3 ) 

for any a > 0. Setting c$ (p, v) := [\/2n N\ N^ 2a " l ^ G exp(2iV + (/x + i/JiV 1 / 3 )) -1 , 
we therefore obtain 

lim ct ( " ) (2^ + /^- 1 / 6 ,2v / A^ + WV~ 1/6 ) 

In (mat, ^jv) - 9n(^n) 9n{vn) 



N^oo I / 



lim 



Cp{n,v) (InHi^N^n) - 9n(Vn) 9n{vn) 



limAT^oo c$ (fi, v) fffl (fi N , v N ) 
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where we have used Proposition 12,11 as well as the assumptions j( a '(/i, //) > 0, 
/( Q ) (z/, v) > 0. This completes the proof of Proposition 12.31 □ 



Proof of Proposition \2.J\ First of all, note that the definition (|2.2p may be extended 
to the case a = and that 1^ (x, x) = Ai(x) 2 for any xGlby Lemma 13. 11 Thus, 
(x, x) > for any i£l, with strict inequality for x > (since it is well-known 
that the Airy function does not have any zeroes on the positive half-axis). More- 
over, note that for any a > 0, 

V / 5x V4tt 3 / 2 y_ 00 (1 - m) a +(V2) y 



1 f +oc exp (^(l-iu) 3 - x(l-m)) 

' - — (in 



4^ 3 / 2 7-oo (1 - iu) a ~0-/ 2 ) 

for any Since linx^oo I^ a \x, x) = 0, this implies that for any a > 0, 

I^(x,x)= / 



for any i£l, Proposition l2.4l now follows by a straightforward induction on a. □ 
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